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SUMMARY 


Accurate calculations to evaluate the perforrnance of the 
stoichiometric carbon-oxygen propellant system have been car- 
ried out for nozzle flow with and without chemical reactions and 
with and without vibrational adjustment. The calculations show 
that, for frozen chemical flow, a lag of vibrational energy states 
at chamber conditions nearly doubles the reduction in Isp, as 
compared with flow in which complete vibrational equilibrium is 
rnaintained. On the other hand, lags in vibrational adjustment 
have practically no effect on the theoretical performance of hot 
propellant systems if chemical equilibrium is maintained during 
nozzle flow. The preceding conclusions are in agreement with 
the results on other propellant systems obtained previously by 


use of an approximate evaluation procedure. 
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I DITRONUCTION 

It is well-!mown that the overall performance of rocket pro- 
pellant systems is dependent on the extent to which physico-chemical 
changes occur during nozzle flow. Thus it ia customary to report 
specific impulse data for "equilibrium flow" (the chemical reactions 
are fast enough to maintain thermodynamic equilibrium at all times) 
and for "frozen flow" (the chemical composition remains unchanged 
at the concentrations corresponding to thermodynamic equilibrium 
in the combustion chamber). Practical experience has shown that 
this method of estimating reasonable upper and lower bounds for per- 
formance evaluation is generally in accord with experimental resulte. * 

Implicit in the customary propellant evaluation procedures is 
the assumption that complete thermodynamic equilibriutn is main- 
tained during flow with respect to excitation of internal degrees of 
freedom. [tis apparent, however, that sizeable departures from 
equilibrium internal excitation may occur, for example, for the vi- 
brational degrees of freedom, particularly for non-hydrogen contain- 
ing propellant systems in small motors for which the vibrational re- 


laxation times may not be short compared to the residence time. (1, 2) 


2eoe @ w= 


«The usual simplifying assumptions in propellant evaluation, which 
do not involve consideration of rate processes, are the following: 
Thermodynamic equilibrium is reached in the combustion chamber 
after adiabatic reaction; expansion of the combustion products through 
the Laval nozzle is adiabatic; the products of combustion behave as 
ideal gases; the adiabatic expansion may be considered to involve 
one-dimensional flow (parallel to the nozzle axis) of non-viscous ideal 
gases; the velocity of the gases at the nozzle entrance position is neg- 
ligibly small compared to that at the exit position. ior further details 
concerning customary evaluation techniques see, for example, ‘'Quan- 
titative Lvaluation of Rocket }ropellants'" by 5. S. Penner, American 
Journal of Physics, Vol. 20, 1952, pp. 26-51. 


.:pproxiinate estimnates of the effect of vibrational! excitation on 
performance of rockets have been carried out previously. (3) 

It is the purpose of the present analysis to present the first 
accurate evaluation of the probable combined effects of chemical 
and vibrational lage during nozzle flow. Ferforrnance calculations 
have been carried out for the following special cases: 


(a) Complete thermodynarnic equilibrium is maintained 
(“equilibrium flow'’). 


(b) Complete thermodynamic equilibrium is maintained 
with respect to internal energy states but no chemical 
reactions occur during nozzle flow ("frozen flow"). 


{c) Complete thermedynamic equilibrium is maintained 
except that no adjustments of vibrational degrees of 
freedom occur during flow (chemically-~equilibrium 
and vibrationally-frozen flow). 


{d) Complete thermodynamic equilibrium is maintained 
except that no chemical reactions and no adjustment 
of the vibrational degrees of freedom occurs during 
flow (chemically- and vibrationally-frozen flow). 





LD. BASIC JATIONS 


fhe requisite basic equations for cases (a) to (d) have been 


a 
piven elsewhere. (2) I or the present perpuses it will be convenient 


to express the cherrical composition by the mole fractions for the 


various chemical species. In terms of mole fractions it is readily 


Shown that equations (7), (S), (11) and (12) of reference (4) become, 


respectively: 


(2) For equilibrium flow 


n n 
La Oe ee 
Bote * ait. a = Nie - oke am. / ey Nie "s e (ie) 
(b) For frozen flow 
VY ec ee ; 
hd he Fat “he rie 7 Hone | (1b) 


(c) For che:nically-equilibrium and vibrationally-frozen flow 


NO 


L = 5 ee oc * n vib .,vib 
= Moue s aHe® = AHCe +(M_/M,) 2, NP vil ee a | de) 


(i) kor chemically~ and vibrationally-frozen flow 


po) 


— ae A wee Cc ryvib vib 
ne j 3 Fr a4 a a 
MU, = AHS = ARS o* 2 Wee (F live Me ) (1d) 


a4 


In equations (la) to (1d), U, represents the linear flow velocity 


° ° s or) '@) s 8 e 
at the nozzle exit position, i and Nye are the equilibrium col. frac- 


tions of species k evaluated at the chamber exit and nozzle exit posi- 


tions, respectively; Uh ke and Mo . equal the total enthalpy (sensible 


“Tor the relatively simple case of flow with frozen vibrational energy 
states, it is not necessary to use the analysis given in reference (4) 
since equations (la) to (id) can be written down immediately. 





plus chemical enthalpy) per Pole of @pecies 4k Uf complete internal 
equilibriuin exists at the nozzle entrance and exit positions, reapec- 
tively; the superscript vib signifies that only the vibrational enthulpy 
is involved; the subscript o denotes complete energy equilibriun. 
exists for all constituents; the superscript o denotes complete chom- 
ical ecuilibrium exists for all constituents. 

I-rom equations (1a) to (1d) it is apparent that the evaluation 
of u, and hence of Isp = u/s requires the calculation of equilibrium 


and N? as well as of total and vibrational en- 


re) 
ke ke 


thalpies. Since the exit temperature is determined by the assump- 


mole fractions N 


tion that the flow process is isentropic, it is obvious that the exit 


re) I qvil 


temperature Th and hence the nurnerical values of Nee’ is and 


Fike? will depend upon the assumed type of flow process. 





Ili, Url. CALCULATIONS 


“or the stoichiometric solid carbon-paceous oxygen pro- 
pellant system [3 mole of Cf{c) reacting with 1 rnole of O,(g)} the 
equilibrium mole fractions at any temperature may be determined 


conveniently by writing the overall chernical reaction in the form 


C + O,~*(1-x) CO, +x CO + [(x-y)/2] 0, + yO (2) 


we arrive at equation (2) by assuring that the first chem- 
ical reaction taking place in the combustion chamber is C+O0,-" CO, 
which is followed by a dissociation of part of the carbon dioxide to 
carbon monoxide and oxygen according to the relation xCO,— xCu 
+ (x /2) O>- This decomposition is then followed by a dissociation 
of part of the molecular oxygen to atomic oxygen according to tire 
relation (y/2) Oo, yO. In equation (2) x is the number of moles 
of carbon dioxide dissociated and y/2 is the number of moles of 
molecular oxygen dissociated. 
From eguation (2) we have 
ny, =z=(2+x+y)/2 (3) 
where n, = total number of moles of reaction products. 
For the chemical reactions under discussion two equilibrium 
constants are obtained ‘ 
aj 
KS Po (Po, 


i 
S32 * Poo (Po MPco,) = Kg /iny (5) 





f 


om ea 


where ov i. and | are in the notation of data cornpiled frotn 
4 


10 
reference (5)* and P). is the partial pressure of species hk. 
For an ideal gas 


Py * mn. Py (ny (6) 


where ny is the number of moles of specics k and Ps is the total 
pressure. 


By use of equation (6) in equation (4) we obtain 


¥ 
Ky = y [ (x-y) /2] cA 2p, M2 +a + yy" (7) 


and, similarly, equation (5) becomes 
1 l 
Ky = « {Lov ix) | 2p, M2+x+y)] 2 7) 


i quations (7) and (5) form a system of two equations in two 
unknowns, x and y. By solving the equations for the expression 
l- 


[ p(x-y) /(2 +x + y)] 2 and equating the results it is readily shown 


that 
x” (p,-KS + K,K,) + (315 Z K,K)-2k5 e 0, (9) 
and 
2 : 1 
y = Kx*/ [ s,-2K,) “ 2u\ (10) 


From equations (9) and (10), the quantities x and y, and the 
number of moles of each of the cornbustion products, can be deter- 
mined” provided Py and T a the adiabatic bak temperature, are 
known. 


* Class notes for JP 200, California Institute of Technology, 19251. 
*%* The effect of lag in vibrational energy states on the equilibrium 
constants has been neglected. It will be shown in Section WIC that 
the changes in equilibrium constant will not materially affect the 
numericai values obtained in the present calculations. 
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In the calculations for finding 1, the following relations are 


used 
se re) oO re) 
Q., 2 -(n. 4H +n Sti +n .At, ) (11) 
AV CoO fe6 CO, “Gas O fy 
= t - 
her o Me LMoxe ” Hoxzge. 16°) my 
and 
QQ, 72. (8) 


where —™ is the heat reieased by the chemical reaction, Ce is 


Oo 


fk ia the stand- 


the heat absorbed by the combustion products, and AH 
ard heat of formation of species k. 

With equations (9) to (13) it is found by trial and error that 
the adiabatic flame temperature is 37 50°R for P, = 20 atmos. (6) 

The mole fractions as found by use of equations (9) and {10) 
and the average molecular weight M, as calculated from the reiation 


M = 2 NM, (14) 


c 
are piven at 3750° and P, = 20 atmos. and at various temperatures 


for P, = latmos. in Table I. Here No and M,. are the mole fraction 


and molecular weight, respectively of especies k. 


TABLE 1, EQUILIBRIUM MOLE FRACTIONS No AND AVER/.GE 
MOLECULAR WEIGHTS M AS A FUNCTION OF PRES- 
SURE AND TIIMPDRATURE FOR THE STOICHIOMET- 
RIC CARBON-OXYGEN PROPSLLANT SYSTEM, 


ras 
Oe f 
, a 
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A. Ugquilibrium Flow 
The evaluation of the performance for equilibrium flow re- 
quires the introduction of an additional assumption, namely, isen- 
tropic flow during the expansion process through the nozzle. This 
condition can be expressed in terms of the component entropies by 


the relation 


[Oo se 
= Neel 
Ix 


0 


xelS 


Oo 
Glee -Rinp,, -RinN,.) (15) 


O yeh / 
oke 7 inp, -RinNy, ee, /Mi.) 2. N 


where Soke and Soke denote the total molay enthalpies at the cham-~ 
ber and exit conditions, respectively and R is the molar gas con- 
6tant. 

The exit temperature, T e can now be found by computing 
the right-hand side of equation (15) for several temperatures and 
comparing the results with the value of the left-hand side. Dy using 
the values of M and Ny as given in Table I and the entropies tabu- 
lated in reference (5), the left-hand side of equation (15) is found to 
be 67.65 cal. /"Kemole for T, = 3750°K and p, = 20 atmos. and the 
right-hand side is found to be 66.43 cal. /°K-mole and 68.16 cal. /°- 
mole at 2900°K and 3000°K respectively for P. = latmos. Sy linear 
interpolation it is found that T e = 2970°K. 

By use of equation (la) it is found that AH?® = 22.600 keal/ 
mole and 17. 437 keal/mole at 2900°K and 3000°K, respectively. 
Linear interpolation yields the result ALI” = 19.04 keal/mole at 
T. = 2970° 2. whence it follows that u, = 2,200 meters /sec. From 
the relation 


¢ 
Lap = u,/8 (16) 





LL). 


the specific imipulge is seen te have 2 value “if $27 seconds. 


>» Frozen Flow 


l or frozen flow equation (15) eimiplifies to 


a - Sone) = Rinip, /p.) (17) 


‘as © allie 


. ni = @) j a nee | ; ; 
since Ni. Nee and ! i. =M.. By proceeding as in Section HIA 
we find 
Te = 2110°K; 


from: equation (1b) 


Cc 


AW” 217.12 keal/mole 
oe 


ig obtained; finally ues 2105 meters /sec. and Isp = 215 sec. 


C. Chemically~Eguilibrium, Vibrationally-F rozen I low 


For chemically-equilibrium but vibrationaily-frozen flow, 


the isentropic relation which is used to determine T ‘ becomes 


FNC. (Sop. - Rinp,-RinN? |) = (Mi, /) 
k= 


j—e 


138} 
= ie vib .vib 
ae ke «ih “3 ce inp, “Rinne. (Sieg ~ "ke | 


ke] 
For diatomic molecules the vibrational contribution to tne 


(7) 


entropy is, in goot! approximation, 


bs -1 
ee eee ee | 


Here u = hcw* /kT, h represents Planck's constant, c is the velocity 
of light, k the Boltzmann constant, and w * the wave nurnber corres- 
poirling to an energy transition from the ground vibrational energy 
level to the first excited vibrational state for the non-rotating diatom- 


ic molecule. For the linear triatornic molecule COs, the total vibra- 





wo | () oe 


(7) 


tional entropy is gi en by the approximate relation 
A , 
vib ae 
Cc — 
=R2Q_ 44, | exp(u,)-1 -In| l-exptu. 20 
Leon l peu;) (20) 


i=] 
where u, corresponds to the four normal vibration frequencies of CO >. 
For the present calculations the numerical values of the characteris - 
tic temperatures 9 = hcw* /k which are listed in Table II were used. 


TABLE NW. CHARACTERISTIC TEMPERATURES FOR O 
CoO 


| 2 


2228 3067 955.6, 955.6, 1988, 3493 


The thermodynamic function lead as found by use of equations (19) 


2? Co, AMD 










and (20), is listed in i:able HI for various temperatures. 


TABLE I, s¥?” (eat. /°K-mole) AS A FUNCTION OF TEMPERA- 


TURE FOR O,, COAND CO, | 

3750°K 3000°K 2900°K. 1600°K 
2228 |.5941|3. 050] 0. 7427 |2. 628 lo. 7683} 2. 5597| 1,393 1.482 
3067 |.8178]2. 442| 2.0223] 2. 027 |1. 0575/1. 9689| 1.9171 6.972 
, Co, | 1988 |. 5301/3. 271] 0. 6627 2.639 |o. 6855 2.7769] 1. 243] 1.677 
460, 1955.6 |. 2548 6. 3185/4. 270 lo. 329514. 2025 0. 597| 3. 0391 


,CO, | 3493 |.9314|2. 198]1.1645]1.774 [1. 2045/1. 7384 2.184] 0. 789 











Com- 
ponent 









lis 
Py use of the appropriate values of m, and M listed in 

Vable I and the appropriate values of se listed in Table HI we 
can now proceed to find Te from equation (13) by a trial and error 
method sirnilar to the method described previously. With the val- 
ues already obtained we find the left-hand side of equation (18) to 
be 67.65 cal/°K-mole at Ts 3750°K, and P. #20 atmos. The 
right-hand side of equation (18) is found to be 67. 557 cal/°K-mole 
and 69. 012 cal/"K-mole at 2900°K and 3000°K, respectively. By 
linear interpolation the result T a 2910°K is obtained. 


For diatomic molecules the vibrational contribution to 


the enthalpy is, in good approximation, 7) 


] 


P =z RTu Lexp(u)-1] ’ (21) 


vi 
H 


lor the linear triatomic molecule CO, the total vibrational en- 


thalpy is given by the approximate relation: 
vib 4 - 
Hy eRTD u, } expla) A] (22) 
iss] | 


Using these relations the values for yyvib {in cal. /mole) lieted in 


Table IV are obtained. 


sy 
TABLE IV. H”*? AS A FUNCTION OF TEMPERATURE FOR 
O,, CO AND Co,. 





LH” (gh voter es 











ae 
\3 has already been noted, tne effect of lag in the vibrational 
energy states on the equilibrium constants has been neglected. In 
order to see just how much of an errvr this assumption introduces 
we will now proceed with the computation of K for T = 3000°K by 
taking the effect of lag in the vibrational energy states into accourt. 
The starting relations for the reaction CO, =? CO = s OO» at the temp~- 


erature T. = 3000°K are 
Al’< =-RYin t (23) 
and 
re ed, ie = i 


from the given chemical re- 
action at 298.16°K) (24) 


where F is the Gibbs free energy as defined by the relation 
FeH- Ts, (25) 


and the superscript e again denotes conditions at the exit temper- 
ature. The enthalpy, Me? and the entropy, Shar are found from 


the relations 


Hie © ene * ke . hic (26) 
vib_ ovib, frye 
a a S ste “(5,07 “ke (27) 


respectively. By use of the enthalpies and entropies tabulated in 


reference (5) and of equations (25) to (27), together with the approp- 


wib 


riate values for S and nn? ub listed in fables HI and IV, the data 


listed in Yable V are obtained. 





> 


= 375 
k( 2 


; j= 3 ae 
12.99 












Sy le ble role mole 


“Te + 9 2 s 
In Table VY, Az ¢ represents the standard free energy of formation. 








Eratn column (5) of Table V we nave 


yaaa * @ Oo lad 
AT (resulting from the given chemical reaction at 298.16 1) 


= =32.803 ~ (-94. 260) = 4+ 61.452 


and, using equation (24), 


AF® s + G6. 08 keal. /raole. 


Using this value for AF” in eryuation (23) it is readily shown that 


= 


: __ #2 — 
sy =? — = Ta li 1. 0138 
OY 
log K, = - 0.442 at Y= 3000°H. 


In the other hand, for complete eneryy equilibrium, 


-0.4692 at © = 3000°K, 


log oa 


3250°K, 


w 


log v4 = <0. 1o9l at L 





~14- 


anc, ising an appropriate Mnear interpbletion, 


-3 


log Ik, & ~- 0.442 at 1/T = .3314 x 197” of T = 3017°RK. 


“rorn the preceding results it le evident that neglecting the effect 
of lag in the vibrational energy states on the equilibrium constant 
IK, at i = 2000°N hag introduced an error equivelent to 17°K. Rep- 
etition of the preceding calciuJatlons showe that the crror in ry 
will be even smaller. Therefore, we arrive at the conclusion£ 
that the error in the mole fractizns at Te = 2910" is small enough 
So that the resultant Te for this case ie within 15° of the value whic” 
would be ovtained uf the effect of lug in the vibrational energy states 
on the equilibrium constants bad »een taken into account. 

Using the calculations given for case (2) it is found that 
AHO’ = +22. 26 keal/mole at 1, = 2910°K from which it follows, by 


vid 


using appropriate values of N? and M from Table I and of I from 


Table IV in equation (Ic) that 5 Mul = AHO° = 4 18.72 keal/mole, 


Cc 


u, 2,203 meters/sec, and a. = 225 geconds. 
& 


S. Constant-Composition and Vibrationally-lrozen F low 


‘or constant-composition and vibretionally-frozen flow tie 
isentropic relation, which is used to deter: ine Tee becomes 


n 
O vib _ ib 
a ne) S.., - iar - (S a Si >] = Rin(pc/pe). (28) 


With the values of Me ¢ Usted in Table I, the appropriate values of 
ao given in Table IN, and the entropies tabulated in reference ( ), 


it is found from equation (23), by use of a trial and error :nethoc, 





=18- 
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ace 
the appropriate values of Fi Y* given in Table IV, and th 


thet T, = 1600°%. Using the values pf Ul tabulated 
' 
tabulated in reference (5), in equat 


Ae = 15.144 kealj« ole, 


J = : . —_we 
. - ; : ze +4 ¥ . S| 2 tl, = 
“, — Le 981 aa a U ww & i j “ Ne . : al ad ; Ba’ : 
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IV. DISCUSSION OF RESULTS 


for convenience the results of the present analysis arc 
summarized in Table VI. Reference to Table VI shows that the 
effect of vibrationally-frozen flow is of no consequence when chem- 
ical equilibrium is maintained. In this case lag in vibrational ad- 


justment lowers the performance by roughly one percent. 


ABLE VI. SUMMARY OF RESULTS 
Chemically-equilibrium and 29°70 


energy-eqguilibriuin flow 





















Chemically-frozen and 
energy-equilibriuin flow 










Chemically-equilibriurn anc 
vibrationally-frozen flow 











Chemically-~-frozen and 
vibrationally-frozen flow 


On the other hand, the effect of vibrationally-frozen flow 
for frozen chemical flow is such as to decrease Isp approximately 
6°/o. In this connection it is of interest to note that lags in vibra- 
tional adjustment, for a given propellant system, become more 
likely as the nozzle size is reduced because a reduction in nozzle 
eize means increased cooling rates. In particular, it has been 
postulated that frozen-vibrational flow occurs in the emall nozzles 


used in the pneumatic method for measuring temperatures in 





(3) 


combustion chambers. ‘~ 
The calculations on the carbon-oxygen propellant system 


(3) 


have been repeated by using an approximate procedure, which 
permits ready estimates of the effect of vibrational lag on per- 
formance for propellant systems for which the usual theoretical 
performance data are available. The results of these calculations 
are surmmarized in Table VII, together with data for other pro- 
pellant systems which were obtained by utilizing the concept of a 
constant average vibrational heat capacity during flow. (3) Refer-~- 
ence to Sable VII shows that lag in vibrational energy states pro- 
duces a similar effect on performance for widely different propel- 
lant systems. A simple physical explanation for this observation 
is obtained by noting that frozen vibrational flow reduces the ef- 
fective heat capacity during expansion by nearly the same extent 
for chemically-frozen and for chemically-equilibrium flow. On 
the other hand, the effective heat capacity is much larger for flow 
in which chemical changes occur provided the propellant system 

is suificiently hot. Hence the effect of vibrational lags will always 
be more pronounced for flow without chemical change than for 
chemically-equilibriurn flow. Furthermore, the disparity will be 
largest for the hottest propellant systems. The preceding remarks 
are obviously in accord with the data given in Tible Wi which show 
a large effect of vibrational lag on performance for chemically- 
equilibrium flow only for the reletively cool RF NA-Aniiine propel- 


lant systern. 





TABLE. II. CF WE IN VIGRALTIONA LeGY siAt 
. ry ANC , VAAIOU ROe"7 ah 


Isp (energy eqailibriwa flow). 


Yropeilant Systeisr. Isp (vibrationally frozen flow) for 


quilibrium § low rozen Pf low 


Carbon-Oxygen* 1. 064 





* Accurate calculations. 
%* /spoproximate calculations. 


“#%*%* Results given in reference (3). 
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